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Abstract 

^^ , An important problem in gas and fluid dynamics is to understand the behavior of 

^^ ' vacuum states, namely the behavior of the system in the presence of vacuum. In par- 

ticular, physical vacuum, in which the boundary moves with a nontrivial finite normal 
P I ' acceleration, naturally arises in the study of the motion of gaseous stars or shallow water, 

•^r , Despite its importance, there are only few mathematical results available near vacuum. 

r ' The main difficulty lies in the fact that the physical systems become degenerate along the 

'^ I vacuum boundary. In this paper, we establish the local-in-time well-posedness of three- 

2 ' dimensional compressible Euler equations for polytropic gases with physical vacuum by 

considering the problem as a free boundary problem. 



^ 1 Introduction 

TJ" ■ Compressible Euler equations of isentropic, ideal gas dynamics take the form in Eulerian 

''^ , coordinates as follows: 

in 

o; 5t/9 + div(H = , 

dt{pu) + div(/9n (^u) + gradp = . 

Here p , u and p denote respectively the density, velocity, and pressure of gas. In considering 
^ ' the polytropic gases, the constitutive relation, which is also called the equation of state, is 

^ . given by 

p = Kp^ (2) 

where K is an entropy constant and 7 > 1 is the adiabatic gas exponent. 

Compressible Euler system ([1]) is a fundamental example of a system of hyperbolic con- 
servation laws. The first and second equations express respectively conservation of mass and 
momentum. It is well-known [15| that the system ([1]) is strictly hyperbolic if the density is 
bounded below from zero: p > 0. When the initial density function contains a vacuum, the 
vacuum boundary F is defined as 

r = c/{(t, x) : p{t, x) > 0} n d{(t, x) : p{t, x) = 0} 
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where d denotes the closure. It is convenient to introduce the sound speed c of the Euler 
equations (P) 
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Kjp 2 for poly tropic gases) 

We recall that for one-dimensional flows, u±c are the characteristic speeds of the system. If 
r is nonempty, the Euler system fails to be strictly hyperbolic along F, namely degenerate 
hyperbolic. 

A vacuum boundary T is called physical if the normal acceleration near the boundary is 
bounded from below and above: 

-oo<^<0 (3) 

in a small neighborhood of the boundary, where n is the outward unit normal to T. In other 
words, the pressure or the enthalpy (= c^/(7 — 1)) accelerates the boundary in the normal 
direction. This physical vacuum can be realized by some self-similar solutions and stationary 
solutions for different physical systems such as Euler equations with damping, Euler-Poisson 
or Navier-Stokes-Poisson systems for gaseous stars [STJ [28l HBJ [63] . 

Recently, the physical vacuum has gotten a great deal of attention among mathemat- 
ical community (see [Ij and the review paper [30]). Despite its physical importance, the 
local existence theory of classical solutions featuring the physical vacuum boundary even for 
one-dimensional flows was only established recently. This is because if the physical vacuum 
boundary condition ([3]) is assumed, the classical theory of hyperbolic systems can not be 
directly applied ^46^i63j: the characteristic speeds u±c become singular with infinite spatial 
derivatives near the vacuum boundary and this singularity creates an analytical difficulty 
in standard Sobolev spaces. Local existence for the physical vacuum states was given by 
the authors f29]. In f29], the authors consider the one-dimensional Euler equations in mass 
Lagrangian coordinates. Existence was proved using a new structure lying upon the phys- 
ical vacuum in the framework of free boundary problems. The nonlinear energy spaces in 
|29j are designed to guarantee such minimal regularity that the physical vacuum ^ is re- 
alized. Coutand and Shkoller fT2] constructed i^^-type solutions with moving boundary in 
Lagrangian coordinates based on Hardy inequalities and degenerate parabolic regularization. 

For multi-dimensional flows, parallel to a lot of activity and important progress in free 
surface boundary problems, in particular regarding incompressible fluids [21 [3l [U dU EOj [351 
[58l \6T\ I64j . more recently, there are some works trying to prove the local well-posedness of 
physical vacuum in 3D. Coutand, Lindblad and Shkoller p^ established a priori estimates 
based on time differentiated energy estimates and elliptic estimates for normal derivatives 
for 7 = 2. As noted in [14], in order to carry out the additional elliptic estimates, sufficient 
smoothness of solutions was assumed and its justification by their energy should require ad- 
ditional work. The purpose of this article is to provide a new analysis of physical vacuum 
based on a hyperbolic type of weighted energy estimates of tangential and normal derivatives 
and to establish the local well-posedness of the three-dimensional Euler system with free 
surface boundary which moves with nonzero finite acceleration towards normal to the bound- 
ary. Independently of this work, Coutand and Shkoller [13j extended their one-dimensional 
methodology and combined it with the a priori estimates given in [H] to construct smooth 
solutions of the three-dimensional Euler equations. We were also informed that Lindblad 
|39j has a similar result using the linearized compressible Euler equations with a Nash-Moser 
iteration. The methods are very different. 



1.1 Existence theories of compressible Euler flows 

Before we formulate our problem, we briefly review some existence theories of compressible 
flows with vacuum states from various aspects. We will not attempt to address exhaustive 
references in this paper. In the absence of vacuum, namely when the system is strictly 
hyperbolic everywhere, one can use the theory of symmetric hyperbolic systems developed 
by Friedrichs-Lax-Kato [19l [331 ES] to construct smooth solutions; for instance, see Majda 
[38] . The breakdown of classical solutions was demonstrated by Sideris [59] . 

When the initial datum is compactly supported, there are at least three ways of looking 
at the problem. The first consists in solving the Euler equations in the whole space and 
requiring that the system ([T]) holds in the sense of distribution for all x € M and t € [0, T]. 
This is in particular the strategy used to construct global weak solutions (see for instance 
Diperna p6] and [10^ I40j ) . The second way consists in symmetrizing the system first and 
then solving it using the theory of symmetric hyperbolic system. Again the symmetrized 
form has to be solved in the whole space. The third way is to require the Euler equations to 
hold on the set {(t, x) : p{t, x) > 0} and write an equation for T. Here, the vacuum boundary 
r is part of the unknown: this is a free boundary problem and in this case, an appropriate 
boundary condition at vacuum is necessary. 

In the first and second ways, there is no need of knowing exactly the position of the 

vacuum boundary. DiPerna used the theory of compensated compactness to pass to the limit 

weakly in a parabolic approximation of the system and recovered a weak solution of the 

Euler system (see also [30] where a kinetic formulation of the system was also used). Makino, 

Ukai and Kawashima [52] wrote the system in a symmetric hyperbolic form which allows the 

density to vanish. The system they get is not equivalent to the Euler equations when the 

density vanishes. This special symmetrization was also used for the Euler-Poisson system. 

This formulation was also used by Chemin ^ to prove the local existence of regular solutions 

in the sense that c,u (z C([0, T);H"^{W^)) for some rn, > 1 + d/2 and d is the space dimension 

(see also Serre |56j and Grassin [23], for some global existence result of classical solutions 

under some special conditions on the initial data, by extracting a dispersive effect after 

some invariant transformation). However, it was noted in [491 I50j that the requirement that 

c is continuously differentiable excludes many interesting solutions such as the stationary 

1 
solutions of the Euler-Poisson system which have a behavior of the type p ^ \x — xq\-i-'^ , 

namely c^ ~ |x — xo| near the vacuum boundary. Indeed, Nishida in [55] suggested to 

consider a free boundary problem which includes this kind of singularity caused by vacuum, 

not shock wave singularity. 

For the third way: the free boundary problem, we divide into a few cases according 
to the initial behavior of the sound speed c. For simplicity, let the origin be the initial 
vacuum contact point (xq = 0). And let c ~ \x\^. When /i > 1, namely initial contact to 
vacuum is sufficiently smooth, Liu and Yang |45j constructed the local-in-time solutions to 
one-dimensional Euler equations with damping by using the energy method based on the 
adaptation of the theory of symmetric hyperbolic system. They also prove that (? can not 
be smooth across F after a finite time. We note that in these regimes there is no acceleration 
along the vacuum boundary. For < h < 1, the initial contact to vacuum is only Holder 
continuous. In particular, the corresponding behavior to /i = 1/2 is the case of physical 
vacuum [12l [131 129] . For < /i < 1/2 and 1/2 < /i < 1, its boundary behavior is believed 
to be ill-posed; indeed, we conjecture that it should instantaneously change into the physical 
vacuum. However, there is no mathematical justification available so far. 

The case /i = is when there is no continuous initial contact of the density with vacuum. 



It can be considered as either Cauchy problem or free boundary problem. An example of 
Cauchy problem when h = is the Riemann problem for genuinely discontinuous initial 
datum (see for instance [71 [23]). An example of a free boundary problem when h = is the 
work by Lindblad [38] where the density is positive at the vacuum boundary. 

Having the local existence theory of vacuum states, the next important question is whether 
such a local solution exists globally in time or how it breaks down. The study of vacuum free 
boundary automatically excludes the breakdown of solutions caused by vacuum, which is one 
possible scenario of the breakdown for positive solutions to compressible Euler equations ([l|). 
It was shown in [l3] that the shock waves vanish at the vacuum and the singular behavior 
is similar to the behavior of the centered rarefaction waves corresponding to the case when c 
is regular |46j . which indicates that vacuum has a regularizing effect. Therefore it would be 
very interesting to investigate the long time behavior of vacuum states. 

When there is damping, based on self-similar behavior, Liu conjectured [42] that time 
asymptotically, solutions to Euler equations with damping should behave like the ones to 
the porus media equation, where the canonical boundary is characterized by the physical 
vacuum condition ([3|). This conjecture was established by Huang, Marcati and Pan [26] in 
the framework of the entropy solution where the method of compensated compactness yields 
a global weak solution in L°°. But in their work, there is no way of tracking the vacuum 
boundary. It would be interesting to investigate the asymptotic relationship between smooth 
solutions obtained by solving a free boundary problem of Euler equations with damping and 
smooth solutions of the porus media equation. 

1.2 Other interesting vacuum states 

Physical vacuum as well as other vacuum states appear in the theory of other physical systems, 
where we believe the methodology of this article could be applied. Here we briefly present 
some of them; see [30] for more detailed discussion. 

The study of vacuum is important in understanding viscous flows [181 HI] • When vacuum 
appears initially, studying Cauchy problems for compressible Navier-Stokes equations with 
constant viscosity coefficients yields somewhat negative results: for instance, a finite time 
blow-up for nontrivial compactly supported initial density [62] and a failure of continuous 
dependence on initial data [25] . There are some existence theories available with the physical 
vacuum boundary for one-dimensional Navier-Stokes free boundary problems [47], for three- 
dimensional Navier-Stokes- Poisson equations with radial symmetry [28], and for other related 
models [iTl [54]. On the other hand, to resolve the issue of no continuous dependence on 
initial data in [25] for constant viscosity coefficient, a density-dependent viscosity coefficient 
was introduced in [44]. Since then, there has been a lot of studies on global weak solutions 
for various models and stabilization results under gravitation and external forces: see J371I65J 
and the references therein. Despite significant progress over the years, many interesting and 
important questions are still unanswered especially for general multi-dimensional flows. 

The relativistic Euler equations are known to be symmetric hyperbolic away from vacuum 
[51j . A particular interest is compactly supported relativistic flows which for instance can 
be applied to the dynamics of stars in the context of special relativity. Whether one can 
extend the theory of free boundary problems including physical vacuum developed for non- 
relativistic Euler equations to relativistic case is an open problem. It turns out that vacuum 
states also arise in the theory of magnetohydrodynamics (MHD), another interesting system 
of hyperbolic conservation laws arising from electromechanical phenomena [151 122] . Due to 
the interplay between the scalar pressure of the fluid and the anisotropic magnetic stress. 



vacuum states are richer than in hydrodynamics and their rigorous study in the context of 
nonhnear partial differential equations seems to be widely open. 

Lastly, it is interesting to point out the connection between vacuum states of degenerate 
hyperbolic systems and degenerate elliptic and parabolic equations. One of the main difficulty 
of studying the Euler system with a free boundary is that it leads to a degenerate hyperbolic 
equation due to the fact that the density vanishes at the free boundary. For this we need 
techniques coming from degenerate elliptic equations; for instance, see [5]. Also, similar 
problems arise in degenerate parabolic equations for instance, in porous medium equations 
[8], in thin film equations [21j, in the study of polymeric flows |53j . 

In the next section, we formulate the problem, introduce notations and state the main 
result. 

2 Lagrangian Formulation and Main Result 

2.1 Derivation of the system in Lagrangian coordinates 

The boundary moves with a finite normal acceleration under the physical vacuum condition 
([3]) and it is part of the unknown. The vacuum free boundary problem is studied in Lagrangian 
coordinates where the free boundary is fixed. 

For smooth solutions, the Euler equations ([TJ can be written in Eulerian coordinates as : 



p{dt + u-V)u + KVp"' = 0. 



(4) 



Let r]{t, x) be the position of the gas particle x at time t so that 

rjt = u{t, r]{t, x)) for t > and 77(0, x) = x in J7 . (5) 

As in [14] , we define the following Lagrangian quantities: 

v{t,x) = u{t,r]{t,x)) (Lagrangian velocity) 
f{t,x) = p{t,r]{t,x)) (Lagrangian density) 

A = [Dt]]^ (inverse of deformation tensor) 

J = det Dr] ( Jacobian determinant) 

a = J A (transpose of cofactor matrix) 

We use Einstein's summation convention and the notation F,^ to denote the /c* -partial 
derivative of F: d^F. Both expressions will be used throughout the paper. We use i, j, k,l,r,s 
to denote 1, 2, 3. The Euler equations ^ read as follows: 



Since 






Jt = JAlv\j and J(0) = 1 



together with the equation for /, we find that 

where po is given initial density function. Thus, using A^ = J~^a^, dH) reduce to the following: 

povi + Kaf{po^J-^),k = (7) 

along with 

ril = v\ (8) 

Now let w be 

w^Kpl-\ (9) 

Note that -^w is the initial enthalpy. We are interested in smooth initial enthalpy profiles 
satisfying the physical vacuum condition ([3]) : 

w = Oondn{= r(0)) , -u; > in f^ , 

1 (10) 

— (i(x) < w{x) < Cd{x) in Vi , 

where d{x) is the distance function to the boundary. The equation ([7]) now takes the form 

w;%^ + (u;^+"Afj-i/"),fc = 0, (11) 

where 

1 
a 
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We have used the Piola identity ([16]) to get (fTT|) from ([7]). Note that a > and a ^ oo as 
7 — > 1. Since tjI = v"^, the equation ([TT]) reads as an -q equation: 

w''Vtt + iw'^''AU-'h,k = 0. (12) 

Thus the equation (fT^ can be viewed as a degenerate nonlinear acoustic (wave) equation for 
r]. Multiply (J12p by r/J and integrate over fi: the zeroth order energy estimates formally lead 
to the following energy conservation 

T- [ il'W'lv f + aw^+'^J-^/ A dx = denoted hy^ = 0. (13) 

at Jq [2 J at 

Note that due to the vanishing factor u;^+" at the free boundary, all the boundary terms 
from the integration by parts disappeared. This energy conservation is equivalent to the 
more familiar form in Eulerian coordinates: 

dtJn^t)\2'^ ' 7-1] 
which is the conservation of the physical energy. 

Remark 2.1. The degenerate nonlinear acoustic system (J12p is in a sense equivalent to the 
Euler equation @ with physical vacuum. Notice that the initial datum w = Kp^" is a 
parameter in the equation and that one recovers the Lagrangian density by taking f = pqJ~^. 



2.2 Differentiation of A and J 

In order to have sufficient regularity so that the flow map r] is guaranteed to be non-degenerate 
and smooth and in particular the Jacobian determinant J is to be bounded away from zero 
and smooth, we need estimates of sufficiently high order derivatives of r] or v. The minimal 
number of derivatives needed will be determined according to the strength of degeneracy, 
namely vanishing exponent a. Here we present the differentiation of A and J. 

Differentiating the inverse of deformation tensor, since A ■ [Dtj] = I, one obtains 

dtAf^ = -Ay,s Al ; 9,4^ = -A'^dirj',, Af (14) 

Differentiating the Jacobian determinant, one obtains 

dtJ = JA'y,s ; diJ = JAldir,\s (15) 

For the cofactor matrix a = J A, from (I14p and (I15p . one obtains the following Piola identity. 

a',fc = (16) 

2.3 Notation 

For a given vector field F on Vt, we use DF, divF, curlF to denote its full gradient, its 
divergence, and its curl: 

[DF]] = F\, 



divF =F' 



jr 



k 



[CUIIFX = tijkF 

Here eijk is the Levi-Civita symbol: it is 1 if {i,j,k) is an even permutation of (1,2,3), -1 if 
(i, j, k) is an odd permutation of (1, 2, 3), and if any index is repeated. 
We introduce the following Lie derivatives along the fiow map ry: 

[Dr,F]l = AlF\s 
div^F =AIF\, 
[cnil^FX = €ijkA]F\, 

which indeed correspond to Eulerian full gradient, Eulerian divergence, and Eulerian curl 
written in Lagrangian coordinates. When r/ is the identity map ~ for instance, the initial 
state of the flow map is the identity map as in ([5]) - these Lie derivatives are the standard full 
gradient, divergence, and curl. In addition, it is convenient to introduce the anti-symmetric 
curl matrix Curl^F 

[CuAr,F]]=A]F\s-AlP,s 

Note that Curl^F is a matrix version of a vector cur^F and that |Curl^Fp = 2|curl^Fp 
holds. We will use both curl^ and Curl^. 

2.4 Energy and main result 

To state our result, we need to introduce some energy which behaves differently for normal 
and tangential derivatives close to the boundary. We first introduce the tangent vector fields 



and the normal vector field for the initial domain. 

For given smooth initial domain Q, we consider a smooth function (/> : M^ — ?> M such that 
{x,4>{x) = 0} = d^ and {x,(^(x) > 0} = $7, and that |V(?!)(x)| = 1 for all x G $7^ = {x G 
O, d{x) < 6} for some small 6 > 0. In particular (/)(x) = d{x) in Qs- Let T be a finite set of 
(smooth enough) vector fields in il which are tangent to the boundary and M the set formed 
by the vector field C which is normal to the boundary and such that C = V0 for x G fi^. 
We also assume that at each point x in fi the vectors of T and Af span the whole space M^ 
and that the second order differential operator ^«g7-/3*/5 + C*C is strictly elliptic. In the 
simplified case where Q, = T'^ x (0, 1), we can just take d(^ = d^ and T = {di, 82}- 

We define the energies B , C , T> and £ by 

B^{v) =- Yj rf^+°+"|curl^9;f9^^'Pda: 

?Tt|+n=0 
-. . N 

C^iv)=- Yl d^^'^^'^lc^rlrjd^d^vfdx 

\m\+n=0 

— ^ r ^ 

^'^ ./si I I , n 

\m\+n=0 

- ^ r ^ 

£^{ri,v) = - Yl ^"^"l^^^c^P + d^+''+''\D^d^d^r]fdx 



:i7) 



\m\+n=0 

In the above formulae, the sum is made over tti G N^ and d'g' denotes the derivative 
dT'^d^'^ ...d^' where I = CardT. The total energy is defined by 

TE'' = T£'' [i],v) =1'' {i^,v) + B'' {v) (18) 

Since ??(0, x) = x, the total energy for initial data t;(0, x) = uq is given by 

7T^(0) = - i d^+"dx + Y ~ I ^"^"l^^^c^oP + d^+"+"|curia;fa^nop(ix (19) 

J \l III n J \ I 

\m\+n=\j 

We also introduce the function spaces X ' ,Y ' 



X"'^ = {d^F G L^{Vt) : / (i"+"|a^a?F|2(ix < 00 , < [ml + n < 6} 

Jo. 

ya,b _ |^i±^^^ ^ ^2^^) . f d^+<^+''\Dd]^d2Ffdx < 00 , < |m| + n < 6} 

Jn 



(20) 



In 

equipped with the following norms: 
b 



IPI|2 — 



Y /" (i"+"|5^5^F|2(ix ; IIFII^,,., = Y f d^^''+''\DdJ^d^F\'^dx 



|»n|+n=0 |m|+n=0' 

Note that 

-c^ II ||2 2 



We are now ready to state the main result of this article. 

Theorem 2.2. Let a > 6e fixed and N > 2[a\ + 9 6e given. Suppose that Dw E X 
and that the initial energy ()19p is hounded: T£ (0) < oo. Then there exist a time T > 
depending only on T£ (0) and \\Dw\\.^a,N a unique solution {v,r]) G C([0,T];X ' xY ' ) 
to the Euler equation ([8]) and (J12p on the time interval [0, T] satisfying 

T£^i7],v) <2TE^{0) and p-/||oo<l/8. 
In particular, 2/3 < J < 2. 

Theorem 12 . 21 indicates that the minimal nmnber of derivatives needed to capture the phys- 
ical vacuum ([3]) depends on the value of the adiabatic exponent 7. The smaller 7 is, the more 
derivatives are required to overcome stronger degeneracy caused by physical vacuum. Indeed, 
this phenomenon was captured in our one-dimensional result [29] where the analysis is carried 
out in mass Lagrangian coordinates. While there is some similarity between our new analy- 
sis and the previous one-dimensional analysis, besides boundary geometry, there is another 
critical difference regarding the energy functionals: In one-dimensional V, V* framework, the 
energy space is very nonlinear in that it is not clear at all to deduce some equivalence to the 
standard weighted Sobolev spaces and moreover, the number of V, V* to define the energy 
space is rigid. On the other hand, in the current analysis, the energy functionals given in (j24p 
are indeed equivalent to the standard linear weighted Sobolev spaces as discussed in Section 
13.11 and also the higher regularity can be readily established. 

The method of the proof is based on a hyperbolic type of new energy estimates which 
consist in the instant energy estimates and the curl estimates. As soon as we linearize the 
Euler system, we start to see geometric structures: the full gradient, divergence and curl of 
flow map r] in the equations. The new key is to extract right algebraic weighted structure for 
the linearization in the normal direction such that we can directly estimate normal derivatives 
via the energy estimates: each time we take normal derivative, we obtain more singular 
(degenerate) weight while the main structure of the equation remains the same for tangential 
derivatives. The effect of taking one normal derivative is worth gaining a half derivative. The 
new energy estimates provide a unified, systematic way of treating all the spatial derivatives. 
In the instant energy estimates, the curl part comes with undesirable negative sign. This can 
be absorbed by adding the curl estimates which will be obtained separately from the curl 
equation. For the construction of solutions, we implement the linear approximate schemes for 
G = [X]«£7-(5/3)*9/3 — w~°'d(^w^^'^d(^ + X]r] with A big enough and we use a relaxed curl of G. 
To build the well-posedness of linear approximate systems, we employ the duality argument 
as done in [29]. Finally, rj is found by solving the above degenerate elliptic equation. 

In the sequel we will give a detailed proof in the simplified case = T^ x (0, 1) and then 
indicate the changes to be done in the general case. The rest of the paper is organized as 
follows. In section [3l we rewrite the energies in the simplified case fJ = T^ x (0, 1) and prove 
some Hardy type inequalities. In Section HI we establish the a priori energy estimates. In 
Section [SJ an approximate scheme is provided and Theorem 12.21 is proven. In Section [6l the 
general domain case is discussed. In Section [71 we conclude the article with a few remarks. 



3 A simplified domain 

For simplicity of the presentation, we will present the full proof in the case when the initial 
domain is taken as 

O = T^ X (0, 1) 

where T^ is a two-dimensional period box in xi,X2. We will then present the main changes 
to be done in the general case. The initial boundary is given as 

r(0) = {x3 = 0} U {x3 = 1} as the reference vacuum boundary. 

The moving vacuum boundary is given by T{t) = 77(t)(r(0)). 

We use Latin letters i,j,k, ... to denote 1, 2, 3 and that we use Greek letters f3,K,a to 
denote 1 , 2 only. Recall the weight w - initial enthalpy ([9]) satisfying (llOp . We now rewrite the 
various energy functionals. We use dp^ to denote d^^dl^^ and |7tt.| to denote |7tt,| = mi + m,2. 

We also introduce r] dependent energies. These energies are equivalent to the ones given 
in pT|) with the difference that d is replaced by w and that a factor J~^/" is added. It turns 
out that these versions are more adapted to the energy estimates and have better cancellation 
properties. 

The curl energies are defined by 



N , ^ N 



\m\+n=l \m\+n=l 

TV N 

C^{ri)= V M tyi+^+^J-i/^lcurl^SJ'SJr^pdx^ V C"^-'^ 



(21) 



lr7i|+n=l \m\+n=l 

The divergence energy is defined by 

TV TV 

V^{r])= y — / u;^+"+"J-i/"|div„a^a3"7?pdx= V P™'" 

rn,|+n=l |m|+n=l 

The instant energy is defined by 

TV 

I , 2 jf^ 2 Jq 

\m\+n=l 

TV 

= E+ Yj ^"''" 

|m|+n=l 



f22) 



(23) 



The total energy is defined by 

T£^ = T£^{rj,v)=£^{r,,v) + B^{v) (24) 

Since 7/(0, x) = x, the total energy for initial data ^(0, x) = uq is given by 

TV 

T£^{0) = T£^{x,uo) = E+ V - / u-^+^japjuoP + u;i+°+"|curia;fajnopdx (25) 



\m\+n=l 
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Next we introduce the function spaces X'^'^, y"'*, Z"'^, which wih be used in the con- 
struction of solutions to approximate scheme in Section 15.11 associated to our various energy 
functional (12111 - (123)1: 



X"''' = {u;tF G L'^in) : [ w''+''\d^d^F\'^dx < oo , < lm| + n < 6} 

Jn 

y"'^ = {w^D^F G L'^ift) : I w^+'^+''J-^\D:^d'^d^F\^dx < cxd , < |m| + n < 6} (26) 

Z"'^ = {w^F G L^{Q) : j w^+''+''J-^dfd^F\^dx < oo , < |m| + n < 6} 

Jn 

3.1 Imbedding of a weighted Sobolev space 

For any given a > and given nonnegative integer 6, we define the weighted Sobolev spaces 
F"'''(J7) by 

F"'*(rj) = {dtF G L^{Vl) : I d'^\D^F\^dx < oo , < A; < 6} 

Jn 

with the norm 

b 

||F||^,,6 =^ d''\D''F\^dx 

We denote the standard Sobolev spaces by H^. Then for b > a/2, the weighted spaces H°''^ 
satisfy the following Hary type embedding [34J: 

with 

II 7?ll -< II T?\\ 

As an application of the above embedding of weighted Sobolev spaces, we first obtain the 
embedding of X"'^ into Sobolev spaces for sufficiently smooth w. 



Lemma 3.1. For h > \cP\, 
In particular, for 6 > [a] + 4, 



I Z7II -i II Z?ll 

l-T II „i^ 3 ll-f^ IIX".'' 



:< m\xo 



If ^ = [D?7]~^ is close to the identity - for instance, if \\A — /||oo < 1/8 which implies 
2/3 < J < 2 - then Y"' and Z"' are isomorphic to standard linear weighted Sobolev spaces 
like X"'''. Therefore, Lemma 13.11 dictates the similar embeddings for y"-^ and Z"''': for 
b>\a]+ 1, 

IIDFII b-a-l T^, ||F||va,6 and ||F|| b-a~l r^, ||F||7a,6 
II \l^ 2 ~ II 11-^ II "h 2 "^ II 11^ 

4 A Priori Energy Estimates 

In this section, we will prove the following a priori estimates. 

Proposition 4.1. Let a > be given and let A^ > 2[a] + 9 and Dw G X"'^. Suppose rj and 
V solve 1^ and i^ for t £ [0, T] with TS'^iv, w) < oo and I/Cq < J < Cq for some Cq > 1. 
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We further assume that r] and v enjoy the a priori bound: for any s = 1, 2, and 3, 
lN/2] [7V/2]-l 

Y, k'^/'^SgV,. I + X; \w'^/'(f^dy,s\ <oo (27) 

|p|+g=o \p\+q=0 

Then we obtain the following a priori estimates: 

^^^(r?,^;) <7-(^^(r?,t;),S^(no),Co) (28) 

where T{£^ {ri,v),B^ {uq),Cq) is a continuous function of S'^ {r],v), B^{uo), and Cq. In 
addition, we have the curl energy B bounded 

B^'iv) < S^(uo) + GiS'^iv, v), Co,T) (29) 

where Q{£'^ {r],v),Co,T) is a continuous function of £^ {rj,v), Cq andT. Moreover, the a 
priori assumption (I27p can be justified. 

The proof of Proposition 14. II is based on the following two key lemmas. 

Lemma 4.2. Assume as in Proposition \4 ■ 1\ Then we obtain the following 

|{^^(r?,^;)+P^(r?)-C^(r?)}<J-i(f^(r,,t;),C7o). (30) 

Lemma 4.3. Assume as in Proposition \4- 1\ Then we obtain the following 

j^C''{v)<M£''iv,v),B''iuo),Co). (31) 

Moreover, ()29p is satisfied. 

By adding (pOj) and (J3T]) . the main estimate ([28|) follows. Moreover, the a priori assump- 
tion (I27p can be verified within £ by Lemma l3. II and the standard continuity argument. 

As a preparation of higher order energy estimates, we first compute the spatial derivative 
of ^,^J-i/". 

diiA'lJ-^/") = j-^/^diA'i - ^J-'^^+^^/^A'ldiJ 

= -j-y-A';Atdiv',s -i J-l/"4'=^^,a^77^, 
= -J-'/''A>;At.div\s -^ J-l/"4>1^9zr?^, 
-J-'/^A^[A^div'-,s-A'rdiv\s] 
Thus 

diiA^-'h = - J-'/°^' [Dr,div]l - ^ J"'/"A' diYr,div - J~^/°^^ [Cml^div]'; (32) 

Similarly, the time derivative of A^J~^'°' is given by 

at(4 J-i/°) = -J~^/"A^ [D^v]i - \J-^I''A\ div^t; - r^l^A^, [Curl^t;][ (33) 

We remark that after taking derivatives, namely after linearization, we start to see structures: 
in ([32]) and (f33l) . the first term corresponds to the full gradient, the second to the divergence, 
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the last term to the curl. For n > 1, we write df'lA^J ^'°'] as follows. 

71-1 ,oA\ 

p=i 

It turns out that the curl comes with bad sign in the main energy estimates, and we will 
obtain the estimates of the curl separately from the curl equation ()40p or ()4ip as if it were 
an ordinary differential equation. 



4.1 The proof of Lemma 14.21 : Energy estimates for £^ 

We first claim that when taking d^ of the equation (I12|] , we obtain the following higher order 
equations of different structures depending on n: 

w^+^'d^riit + (u7^+"+" d^[AU~^h),k +u;"+"I°'" = (35) 

where I^'"' are essentially lower-order terms given inductively as follows: 

jO,0 ^ . jO,n ^ 53j0,n-l ^ q^^ . Qn-l^j^^ j-l/a^^^ _^^^ . gn^A^j-1/a^ 

+ (a + n)d3W,k ■ d^-^ [A^J'^^''] forn > 1 
In order to see that, it is convenient to write the equation (j35p in the following form: 

dWtt + ^d^-"' [Afj-'h + (! + « + n)dsw ^g^i^' J" V-] 

+ wd^[AfJ-^/%^ +(1 + a + n)w,^ d^'[AfJ-^/''] + l"'" = 

This form can be readily verified by the induction on n starting from ()12p . The equations 
for general mixed derivatives dT'd^r] read as follows. 

lU^+ngmgn^i^ + («;!+"+" 9^9J[^f J-l/"]),fc +u;"+"I™." = Q (37) 

where J'"'" is given inductively as follows: for |?7i| > 1 

jm,n ^ Q^jm-l,n ^ g^^ . ff^-^ d^iA^-'h^k +(! + « + n)dpW,k -dJ-'d^iAWh (38) 

The degeneracy of the equation lies along the normal direction X3. It is interesting to point out 
that the more normal derivatives we take, the more degeneracy we have, while the structure 
of the equation stays the same for the tangential derivatives djs . 

We now perform the energy estimates of ([37|) for dj^d^v and Dr/dJ^d^r]. By using ^^, 
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we rewrite (I37D as follows: for 1 < |?n,| + n < A^, 



|m|+?i— 1 

\p\+q=l 
+ u;"+"/™'" = 

Multiplying by dj^d^rjl and integrating over Q, we get 
1 d 






|rri|+n— 1 



T 

V ,s 



bl+g=i 



(39) 



/■^;"+"a^aJr/^/'"•"dx = 



We will estimate line by line. The first line in (j39p gives rise to the L^ integral of the full 
gradient Dj^d^d^^r] with the weight y;i+"+" j-i/" and commutators involving dtA^ and dtJ- 

- f w^+^+'^d^d^v^k J-^/"A^t • [D^a^a^r?];dx + — f w^+'^+'^J-^Jt \Dr^dfd'^r]\^dx 

= 4-?™'" + {i) 
dt 

For the second line in ()39p . divnd'a^d^'i] and commutators come out. 

- ^ I w'+-+"d^d^v\k J~^'''A\, ■ div,a^a3>drc + ^ y ^i+-+"j-^ Ji idiv^aj^ag^pdx 
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The third hne can be written as Curl^So^SJ ry plus commutators. 



i>r 

= ~IT /^t;i+"^"J~i/"|Curl^9^9j7?|2dx--^ f w^+'^+^'J-^ Jt\Curl^dJ^d^7]\'^dx 

i>r "^ 

= -Ac™'" + {iii) 

Note that we have used the following anti-symmetrization to obtain the curl structure. For 
any (1,1) tensors E and F 

i,r=l K i>r i<r ) 

= YSK-E[MF:-F:) + Y,EI{FI-Ft) + Y,K{FI-F^ 



i>r i>r i<r 



Y^iEl-E^iFl-Fl 



i>r 



Since \A^^\ and j Jj| are bounded by the energy 8^ given in ([M|) . the commutators (i), (ii), {iii) 



are bounded by a continuous function of £^ : 



N\ 



|(i)| + |(i,)| + |(,i,)|<7-3(^-) 



Next we turn into the fourth and fifth lines in ([39]) which consist of lower order nonlinear 
terms. We present the detail only for the first term and the other can be estimated in the 
same way. We consider two cases: \m\ = and \m\ > 1. For \m\ = 0, the index p does not 
appear and the first term reads as follows: for 1 < q < n — 1 



dl4 ■ (u'i+"+" dllJ-^/'^A^^Al] d^S'',s ),k dx 

+ / 814 ■ (u;1+"+"9|[J-^/°A3^|] d'^-'^if,s ),3 dx 
= (iv) + (v) 
For (iv) where k = 1, 2 

+ / 8^4 ■ u;^+"+"af [J-^/"A;:A'] 9«93 "''r/^, dx 
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Since \w^''^d1vi'^ ,s \ and \w^^ ^^''^dpd'^ rf ,s \ for each q < [n/2] are bounded due to the as- 
sumption (f271) . {iv) is bounded by 

n-l 

g=l 

For [v): 

{v)={l + a + n) I w'dlit ■ u-^+^af [J-i/"^3^|] d^-^rf,, dx 
+ / d'^4 ■ u;i+°+"a^+^ [J-i/"^3^|] ^3 "''r/^, dx 

Since Iw'^'^Sl?]^,^ | and \w^'^~^'''^df^d^' rf\s \ for each q < [n/2] are bounded due to (j27|) . (f) 
is bounded by 

n 
9=1 

Now consider |?7i| > 1. Recall that \m\ + n < N, 1 < \p\ + q < \m\ + n — 1, < \p\ < \m\ 
and < g < n. We first integrate by parts in each x^ direction of the second factor and then 
integrate by parts in x^ (/? = 1 or 2) direction of the first factor. 

I d^d^vl ■ (^'+"+" dpl[J-'/''A';At] dJ-^dtS',s ),k dx 

+ I d'^~^dl4,k ■ «;!+"+" ^5|[J-i/"A^A^] a^+'-^aJ-^'r?'-,, dx 

Note that for A; = 1 or 2, the first factor together with the weight w^°^~^^'''^ is bounded by 
£•771,71 g^j^^ £qj. /j = 3^ ^j^g gj-g^ factor with the weight t(;("+"'+^)/2 is bounded by £;»"-i:"+i. Thus 

from (i27|l again, the above is bounded by 

7TI+71— 1 
^711,71^^711-1,71+1^ Y^ gp+l,q < J-g(^^). 

(*) p+g=i; 

0<p<m,0<g<7i 

Note that since \m\ > 1, n + 1 < N and |?tt,| + n < N, (*) is bounded by the energy £^ . 

For the last line in (j39p . first note that J™'" in (|36p and (|38p consist of terms having at 
most |?7z| + n derivatives of Drj with appropriate weights. Hence as done for the previous 
terms, by using ()27p and the fact of Dw G X"'^, it is easy to see that it is bounded by 

cm,n _j_ cm+l,n—l , c"m—l,n , c'm,n—l 

This completes the proof of Lemma 14.21 
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4.2 The proof of Lemma 14.31 : Energy estimates for C^ 

It is well known that in Eulerian coordinates, the curl of a gradient vector is zero. Here is 
the Lagrangian version. 

Claim 4.4. If uj'' = Alh,r, 

CUrlriLO = 

Proof. 

[curlew]* = eijkA'^j lo^ ,s = eijkA'-{Al,s h,r +Alh,rs ) 

— f:., (—A'^ A^ n"^ A'^h J^A^A'^h \ 



Since both terms are symmetric in j, k, the conclusion follows. D 

The equation for ii in (j6]) is equivalent to 

vi + {l + a)AUwJ~'/"),k = 

and therefore, by the above claim, the curl equation is written as follows: 

curl^at7; = {eijkA'jdtv\s = 0) (40) 

It can be also written 

9i[curl^t;]^ = e^^fc^Mj ^'^s (41) 

In order to estimate C"'", we first integrate (j4ip in time: 

rt 

/o 

For < |m[ + n < A^, we take dg^d"!} first to get 



curlr,dtr] = curluo + eijk / 9i^f v'^,s dr 



5t[curl^5^5Jr/]^ = [cmW^d^uo]' 

+ e,,kdtA'j-dfd^v',s- Y. ^^^k%dlA].dJ-^dt'dtr^\s 

\p\+q>l (42) 

+ E ^^Jk fd''^dldtAldJ-Pd^-%',sdT 
\p\+q>0 •'^ 

The last term in (j42p seems to contain more derivatives than the total energy. However, since 
each term is integrated in time, by the fundamental theorem of Calculus, we can move a time 
derivative of one factor to the other. If |p| + (7 > [N/2\, we write it as 

f S^pdldtA'^ ■ dj-''d!;-%\s dr = ^^dlA'j ■ 5J-^a3"-%^, -dplA'^ ■ dj-^'d'^-'v^, (0) 

- f dldlA^d'^-^dt'dtv^.dT 
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and if \p\ + q < [N/2], 

J 

- fcf^dld^A^-d'^-^d-~''v',sdT 

J 

Note that 9p|^j"(0) = for b| + g > [Af/2] and 9J-PaJ"«??^, (0) = for b| + g < [iV/2] 
since [I?r/](0) = /, and thus there will be no contribution of initial data. Hence, the curl 
equation (Ii2]l can be written as 



dt\<^v.T\dfdlT]\ = [curl5^53"«o]' + 2eijudtA'^ ■ d^d^jf,s 

0<|p|+g<[Af/2] 

- E '^^^^ f ^?I^.^ • d^^'dt'dh^s dr (43) 



|p|+g>[Ar/2] 



E ^^3k I dpId^A^ ■ d-^-'d^'v^s dr 



\p\+q<[N/2] 

Note that each factor in the right-hand side has no more than |?77.| + n + 1 derivatives of 
rj. Hence, in multiplying (|43p by w^^"^"'J^^''^[cuvlr^d]^d^r]Y, integrating over Q, and using 
the a priori assumption (127p as done in the previous energy estimates, we obtain the desired 
energy inequality (j31|) in Lemma 14.31 

The energy bound (p9]) for B^ {v) is obtained directly from P3]) since ^([curl^SS^SgTy]* = 
[cwilrjdV^d^v]^ + eijkdtA^- ■ d^d^r]'^ ,s and each term in the right-hand side of (jl3|) is bounded 
by 8 {rj,v) and Cq. This completes the proof of Lemma 14.31 

5 Existence proof 

5.1 Iteration Scheme 

In this section, we implement the linear approximate scheme and prove that the linear system 
is well-posed in some energy space. 

Let the initial data ??(0, x) = tjq = x and dt'q{0,x) = uq{x) of the Euler equation ([12]) 
be given so that T£ (0) < i? for a constant i? > in (j25]) . We will construct approximate 
solutions r]^{t,x) and dtr]u{t,x) for each nonnegative integer u, by induction satisfying the 
following properties: 

r7^(0, x) = x, dtr]u{0, x) = uo{x) , \\Au - /||oo < C^ for Cu < 1/8 (44) 

as well as 

[N/2] [N/2]-l 

Y, \w'^'dplv/,s\+ Yl k''/'5^5|5t^/,s I < oo (45) 

\p\+q=o \p\+q=o 

where the approximate inverse of deformation tensor and Jacobian determinant are defined 
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by 

A^, = [Dr]^]-^ ■ J^ = det Drj^, . (46) 

The condition ||^j^ — I||oo < C^ for Ci, < 1/8 in ([2]) also guarantees the non-degeneracy of 
the approximate flow map: 

1 1 

2/s < < T < < 2 

' - 1 + Sa + 6C72 + 6C3 - " - i-3C^- QC^ _ 6C3 " 



(47) 



The differentiation of A,y and J,y, which is similar to (I14p and (115p . is given by 
(9^4 '^ — — 4 '^B-.n '■ 4*-<9, 4*^ — — 4 *^<9,t) ^ 4 






OtJu = JuAuldtiiJ.s ; diJu = JyAyldiTjJ, 



Note that Piola identity (fTUj) also holds for A^. a,y'^,k= 0. Having A,y and J^, defined, we 
now introduce the following vector- valued, non- negative and weighted second order linear 
operators Cf,, C^ for a given vector G and a weighted first order linear operator C^ for a 
given anti-symmetric matrix H as follows: 

[CIGY = -{w'+"J,~'/"A,';.A,^,G\s),k 

[J^t Gf = -^(u;i+" J,-V-^,f A,?G^, ),fc (48) 

[/:-^i7]^ = _(^i+"j^~i/"A,^/7;),, 

For the iteration scheme, instead of r], we approximate [—df — c^l — w~"^d^ w^^'^d^ + X]t] = G 
for a sufficiently large positive constant A > and Curl^G = H. For each z^ > 0, consider 
the following approximate system for Gu+i and Hu+i: 



where TZi, and Sy are obtained by replacing -q by rjy in 7^ and 5 given in ([66|) and ([67|) . For 
the convenience of readers, we put the formal derivation of G and H in Section 15.51 Note 
that T^v and S^ consist of lower order terms and moreover, 

117^1,11 YQ,iv-2 < oo and ||5,y|La,]v-2 < oo \i E'^ irjy , dt'qu) < cx). 

And also note that [5;^]^ = — [^j^]]. and thus Ryj^i is anti-symmetric. 

The initial data for Gyj^\ and Ryj^i are inherited from the original system: 



G^+i*(0, X) = -(1 + a)d'iw^j, + \x\ 

5iG,+i(0,x) = {-dl - dl - vj-'^dsw^+^'d-s + X]uo{x) , Hy+iiO,x) = 0. 



(50) 



The unique r/j^+i is to be found by solving a degenerate elliptic equation 



The same goes for dtGy^i. Ay^i and Ji^+i are defined as in (j46l) . We remark that the 
approximate systems ()49p converge to the equations for G and Curl^G respectively in the 
formal limit z^ ^ oo. 
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We define the approximate energy functional Sy+i at the v^^ step: 
^^ 1 /■ 1 /■ -1 



\m\+n=Q 



+ ^ [ w'+^+^^j/^ldiYrt^dJ^d^G.+il^dx 



2a 

1 

- / w 



l+a+n 



— "'"Ilrjr' l|2 _i_"''llr' l|2 -LQllH- l|2 

z z ^^ ^^ 

|m|+n=0 

where Y"' ~ and Z^' ~ denote y°'^-2 g^^^^ ^a,N-2 -^^ (j26p induced by r^i,. Note that by 
Cauchy-Schwartz inequahty, 

<- I w^+'^+''J^^\ar^^dfdlGu+i\^dx + - I w^+''+''J^^\d'^dlH^+i\^dx 
and therefore 

''^li;)/^ l|2 _i_"'^llr' l|2 _1_"'"IIH' l|2 ^r 

2 "^ b ^i' 2 ^i' /f-o\ 

1 ^1 7 *- '' 

<r- wan l|2 ! (" _i_ Mir* l|2 I ' II H" ||2 

— 7T r^t'-^!^+l YQi,JV-2 i" It: i" 7; — I U-^+l L.a,JV-2 i" 7: -"V+1 rya,N-2 

2 b 2a ^v Z ^v 

We now state and prove that the approximate system (|49p is well-posed in the energy 
space generated by E^ under the following induction hypotheses: 

(HPl) f I, < 00 and rj^ and dtrjy satisfy dM]). 

(HP2) The left-hand side of (HSil is bounded by fj^ and £u~i- 



(HP3) £ {rjiy , dtr]u) is bounded by £1, and £^, 



u-l- 



Proposition 5.1 (Well-Posedness of Approximate system and Regularity). Under hypothe- 
ses (HPl) - (HP3), linear system (jl9]) admits a unique solution {dtGyj^i^Gu+i^Hi,j^i) in 
X"' , yjf' ,Zu' . Furthermore, we obtain the following energy bounds: 



Jo 

where TT{£,y-^-i,£i,,£^^i,B^ (uq)) is a continuous function of £^^i, £^, £u-i, B^{uo 

Note that TZ^, € Y^' ~ and S^ € Z"' ~ under (HP3). Proposition 15.11 directb 
from Proposition 15.51 given in Section 15.41 
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In order to complete the induction procedure of approximate schemes, it now remains to 
verify the induction hypotheses (HPl), (HP2), and (HP3) for v + 1. 

By the fohowing Elhptic Regularity lemma, which will be proven in Section [5.31 we obtain 
a unique rj^j^i to the degenerate elliptic equation (|5TI) . 

Lemma 5.2. Let k > be given. For each G € X"'*^, there exists a unique solution u € 
j^a,k+2 ^^ ^^g following degenerate elliptic equation 

[-df -dl- w-'^d3W^+'^d3 + X\u = G 

Moreover, we have 

\\u\\xa,k+2 ;< \\G\\x<^,k ■ 

First (HPl). By Lemma[521 iju+i and dtriu+i constructed in the above satisfy the following 

l|f^t^I^ + l|lx".^ + ll^i^ + l|lv".'V ii \\dtGjyJ^l\\xa,N-2 + ||G;y4.l|La,JV-2 

as well as the initial boundary conditions in (j44p . The boundedness of Ju+i will follow from 
the continuity argument by using the estimate of dtJu+i- Note that from Jacobi's formula 

dtJv+i = dtdetDr]^+i = tr:{adi{Dr]^+i)dtD7]u+i) 

where adj(-) denotes the adjugate of a given matrix. Since dtr]u+i G X"'^ and 7]^+! G Yu' , 
by Lemma [3?H \Dri^^i\ and \dtDri^^i\ are bounded by their X"'^ and Y^' norms and hence 
bounded by £u+i, £u- Thus 

\dtJu+i\ < C'l^+i 

where C^+i depends only on £u+i, £u, and initial data. Since Ju+i{t) = Ju+i{0)+Jq dtJu+idr, 
we get 

1 — Cu+iT < Ju+i < 1 + Cu-^-iT 

Now since Ju+i is bounded away from zero for sufficiently small T, Drj^j^i is invertible and 
thus A^^i is well-defined and the formulation in (I47p is also well-defined. To verify the last 
condition in (j33]) for {u + 1)*^, we use ([^7|) . Since A^^i and dtDrj^^i are bounded by £u+i, 
£y and since ^;,+i(t) = ^i,+i(0) + /^ dtAyj^xdr, 

and thus for sufficiently small T, Cj^+i < 1/8 can be found. 

We move onto (HP2). This directly follows from the embedding as in Lemma l3.H since 
N is sufficiently large iV > 2[a] + 9. 

For (HP3), it suffices to show that for each m-|-n < N , J^ w^'^"^^ J~_^{°'\Drj^^-^dJ^d^rjy^i\'^dx 
is bounded by £u+i and £u- Note that 

Jn 

— / -,,.l+a+"(^ '^+^ \-J-/Q 7-1/"|.4 . «r7 ^A ^fV^ff^r, ,J I'^rJ'r 

Because Jy+i is bounded away from zero and Ayj^i and Drjy are bounded by £u+i and £y, 
(HP3) at {u + l)*'^ step readily follows. 
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5.2 Convergence of the iteration scheme and uniqueness of solution 

In order to prove Theorem 12.21 it now remains to show that r]i, and dtr]i, converge, the hmit 
functions solve Euler equations (l8|) and (fT2|) . and they are unique. 

First, by applying the Gronwall inequality to the energy inequality obtained in Proposition 
15. Ij we can deduce the following Claim: 



Claim 5.3. Suppose that the initial data r]{0,x) = x and dtr]{0,x) = uoix) of Euler equations 
<^T2\i are given such that TS (0) < B for a constant B > 0. Then there exist T > such 
that if for all i < v, £i < 3S/2 for t < T, then £u+i !^ 3B/2 for t < T and in addition, for 
all V, \\Ay - I||oo < 1/8 and 2/3 <Jy<2. 

Thus we get the uniform bound of 8p as well as the uniform upper and lower bounds of 
Jjy. In order to take the limit i^ — > cxd, in view of (jSSp . we define the homogenous energy 
functional £u+i for G^+i and Hyj^i: 

|m|+n=0 -^ -^ (54) 



+ iy" «;!+"+" lapji/.+ipdx} 



Then due to ([53]) and since ||^jy — /||oo < 1/8 and 2/3 < Jy <2, £u+i and £u+i are equivalent: 

1 



l + M, 



■£u+i < £u+i < (1 + M^)£^+i 



where My depends only on Cy and £u. Now by Claim [5^31 £i,+i and M^s have the uniform 
bound over t <T. Therefore, there exists a sequence f^ such that G^^, H^^, and rj^^ converge 
strongly to some G, H, rj. Due to the uniform energy bound, we also conclude that t] and dtr] 
solve dH]) and ([12]) with the desired properties. 

For uniqueness, let {r],v) and (?y, w) be two solutions to ([8]) and (I12p with the same initial 
boundary conditions having the total energy bounds: T£ {i],v), T£ (^, t^) < 2B. Define 
Z{t) by 

Z{t) = ^ fw"\v- vfdx + a I u;i+" J-i/"^ V _ j|2^^ 

|m|+n=l 

From (I37p , the equations for r] — rj read as 



In performing the energy estimates of 77 — r/ as done in the a priori estimates and noting that 
I j-1/q _ J"^/"| < CiB\Drj - Dli\ and |^ - A| < C2B\Dri - Dr/I, one can obtain 

— — < CZ for C depending only on B . 

Since Z(Q) = 0, this immediately yields the uniqueness. 
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5.3 Proof of Lemma 15.21 

Here we give a sketch of proof of Lemma 15. 2i We refer to [5] and [6| for similar regularity 
results of degenerate elliptic problems. We introduce the spaces H = X'^'^ and V = X°'^ 
defined by 

We also define the following scalar product on H, (n, v)h = Jq w°'uv. We define on V the 
following norm 



o a_ o ^^ — ^ II ^^ ||9 II a-\-l 

\v = \\w^u\\^2(^a) 

(3=1 



\u\\v = \\'W2u\\Utr,\ +'^\\w2dpu\\l2^Q'j + \\w 2 (93u||^2(f^) 



Lemma 5.4. V is dense in V . 

The proof is based on an explicit approximation. We define Xn{t) = ^ ^^ ^ — ^ — ^ ~ n 
and Xn{t) = nt - 1 if i < f < ^ and Xn{t) = n - 1 - nt if 1 - f < t < 1 - ^ and Xn{t) = 
for 0<t<-orl < t < 1. For each v £ V, we define Un = Xn{x3)v and f" = pj_ * Un 

the convolution of m„ with the mollifier pj_ = {2n)^ p{w-) where p G C^(M^), p > 0, f p = 1 

2n 

and supp(p) € -8(0, 1). Hence Vn G T>{^) and Vn goes to u in 1/ when n goes to infinity. To 
prove that Vn converges to u in 1/, we have to use the Hardy inequality, namely the fact that 



w"-^] 



n Jn 



This ends the proof of lemma 15.41 

We also define the following bilinear form 

f ^ 

B[u, v] = Xw^uv + ^ dpudp{w"v) + w"'^'^ d^ud^v 

Jn p^^ 

where A is big enough. Notice that B is not symmetric if w depends on the tangential 
variable. If A is chosen big enough then, B satisfies the hypotheses of the Lax-Milgram 
theorem, namely the fact that 

I-B[n, ?;][ < C||n||y||i'||y, and — ||w||y < -B[t',f]. 

Hence, for each bounded linear functional / on V , namely / : V — )• M, there exists a unique 
element u € V such that B[u,v] = f{v) for each v & V. One can try to characterize the 
set V, but we do not need to do it here. For any G & H, f{v) = {v,G)h defines a linear 
functional on V and hence there exists a unique u gV such that B[u,v] = J^w^Gv for all 

veV. 

Using lemma [57il and using a density argument, it is easy to see that u solves 

i-df -dl- -iy-^agw^+^ag + A]n = G (55) 

if and only if B\ui,^v\ = J^w'^Gv for all v G V. Hence, u is the unique solution of ()55p . To 
prove that u is more regular, we use again Lax-Milgram to construct C/g the solution of ()55p 
with the right hand side G/3 = dpG + a^vj-'^d-ilw^+'^d^u] + (1 + ajuj-^agfw^+i^agu]. 
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Indeed, one has just to observe that /^(w) = {v, Gj3)h = 



H 

W ' . . . . W ' 



defines a Hnear functional of V. Then by uniqueness of the solution to ()55p with the right 
hand side Gjs, we deduce that Up = dpu. To be more precise, one has to replace the partial 
derivative with respect to xp by difference quotient and then pass to the limit to deduce 
that Up = dpu € V and that ||9/3n||y < C|[G|[h. Using the equation (|55]) . we also see that 
w~'^d^w^^°^d2,u € H and hence by Hardy inequality, we deduce that 



/ w 

Jn 



'^\w^+''d3u\'^dx <C I u;-°|a3W^+°a3n|2dx < C\\G\\h- (56) 

Jo. 



Using that w~'^ d^w'^'^"' d^u = {l+ajd^wd^u+wdlu and the fact that (|56]l yields that d^u € H, 
we deduce that wd^u G H. 

We now assume that G is more regular, namely that G G X"' for some k > 1. We 
want to prove that u G X°' "'"^. From the previous argument, we know that u G X"'^. 
Proving regularity in the tangential direction is very similar to the previous argument. Indeed, 
differentiating ()55p with respect to /3, we see that dpu solves ()55p with the right hand side 
Gp and due to the extra regularity of G, we see that Gp G H. Hence, we can apply X"'^ 
regularity property to dpu and deduce that it is in X°'^. Of course we can repeat this k 
times and deduce that the tangential derivatives d^u G X"'^. It remains to control the 
normal derivatives. Taking one normal derivative, we observe that Ssm solves 

with the right hand side G3 = d^G + (1 + a)d^wdsu. Moreover, it is clear that d^u G H and 
that G3 G H. We introduce two other Hilbert spaces Ha+i and V^+i where for each 7 > 0, 
H^ and V^ are given by the following norms 

2 

p=l 

In particular, we notice that H = Ha and V = Va- 

Now, we can apply the previous regularity argument with H replaced by Ha+i and X"'^ 
replaced by X""*"^'-^. We can also combine the X3 derivatives with tangential derivatives and 
prove that for all m, < m < /c — 1, we have dT'd^u G X"^"^'^. 

By an induction argument on the number of X3 derivatives, we can finally prove that for 
all m, < m < /c — n, we have dV^d^u G X""*""'^. Hence, we deduce that 

5.4 Duality Argument: solvability for fj49|) 

Proposition 5.5. Let r] and dtf] be given such that £ {ri,dtr]) < 00. For f = f{r]) and 
g = gijj^dtrj) where f G L^(0,T;X"'''), g G L^{0,T; Z"'^), and g is anti- symmetric, there 
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exists a unique solution {dfG, G, H) on (0, T) to the linear system 



w'^dfG + C^G + C^G + C^H = w'^f 



w 



l+a 7- 



J-^/'^dtH = 'u;i+" J-i/^ 



a 



(58) 



G{t = 0) = dtG{t = 0)=H{t = 0) = 
and the solution satisfies 

Moreover, if f e X"'^-^ g G Z°'^-^ i/ien 

/or some constant C that depends only on T£ {r],dtri). 
Let A denote the set 



(59) 



(60) 



A 



ij 



G G°°((0,oo) X n) such that (0, Sic/), V')t= 







Hence, {dtG, G, H) solves (1581) on a time interval (0, T) if and only if for each test function 
(i;^, tp) G A, we have 



/ f G- {w'^djct) + £"(/> + £V) - ^-u^^^" J"^/"/i' : Curiae/. dxdt =11 w'^fcf) dxdt 

[ Aw^+'^J-^/'^H : {dt^ - dtiJ-^^'^)^)dxdt = [ [ -Aw^^T^'''g^dxdt 
J Jo J 



(61) 



We denote 



defined on the core 



V 



w 



"a^G + C^G + C^G + C^H 



-w 



l+aJ-l/aQ^H 



dfG G LiX'^^^dtH G Lt'Z"'^G G Li{V{C^)) 



Hence V can be extended uniquely to a closed operator. Moreover, A C W* , the dual of V, 
and 

VVV V^i+"J-V«[45i^ - 4at(J-V")V; - iCurl^c; 
Therefore, (I6ip holds for each ((/>, ^) G ^ if and only if for each (0, tp) G ^, we have 

cT 



T r /r 

) • V* , 






(62) 



We take (^) G A and denote 



V* 
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The energy estimates for V* with (j52p and (j53p yield that 

'^\\dtm.,o + ^11011^0 + ^iiV'iiko} <c [ 

0<t<T 



111 f^ 

sup {-\\^t^j^^,o + -UWycfi + -ilV'lll^.o} < C / \Mx^,o* + \\^\\l^,o'dt 



where X"'^ and Z'^'^ denote the dual spaces of X"'" and Z"'". Thus the operator V* defines 
a bijection between A and V*{A). Let 5o be its inverse. Hence 

So : V*{A) ^ A given by Q) \ ^ (^ 



and we have 



||(9t(/>, 0,V')llc([O,T];X«.Oxy".OxZ".O) < C"!!! ^ ) !Il1(X«>0'xZ".o*) 



We extend this operator by density to V*(^) and to L^(X"''^ x Z"''' ) by 

Hahn-Banach. We denote this extension by S: 

S : L1(X°'°* X Z"'°*) ^ C7([0,T];X°'° x y"'° x Z"'°) , (^) I ^ ("^ 

Now we want to solve ([58]), namely V(^) = (^i+^J/i/c. ) with G(t = 0) = dtG{t = 0) = 
-fr(t = 0) = 0. This is equivalent to requiring that ([62]) holds for each (0, ip) G ^. Hence, it 
is enough to show that for all (^) G L^(X"''^ x Z'^'^ ), we have 



a 



G\ 




^fy 




■k: 


+ qW 


> ^f?^ - 


- Q*( 


"'"/ ^ „.>,o,-o 


Q* ic fT 



Therefore, it is enough to take (^jj) = S* (_4^^,i+a j-i/a ) where S* is the dual of S, which 
satisfies 

S* : X(0,r;X"'°* X Z°'0*) ^ L~(0,T;X"'° x y"'0 x Z°'0) 

and thus ||((9tG, G, -ff)|jj;^tx)(o,T;X"'Oxy"'OxZ"'0) ^ C'|l(/,5)||^i(xa,0xZ"'")- ^t this stage we do 
not know whether {G,H) is continuous in ()59p . This actually follows from the regularity and 
the density argument. 

The uniqueness of {G,H) will follows from the fact that if / = (7 = in (|58p . then (0,0) 
is the only solution to (f58]l in L°°((0,T;X"'° x Y"'^ x Z"'°). To prove it, consider a solution 
{dtG,G,H) G L°°((0,r;X"'0 x y"'0 x Z"'0) to dMD with / = 5 = 0. We will also make use 
of the duality argument. Indeed, as changing the roles of V and V* and arguing as above, we 
can prove the existence of a solution (0, ip) to the dual problem 

y^i+Q j-i/q[4^^^ _ 45^(j-i/")^ _ iCurl^0] = ^ (63) 

[<P{t = 0) = dt<P{t = 0) = i;{t = 0) = 

for each ($,^) G Li(X"'0* x Z"'0*). Then for each ($,*) G L\X°''°* x Z"'0*), we consider 
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jtp) a solution to (|63|) . Hence, we can write (|62|) with the solution (cp^tp). This yields 

;7(^)(:)"- 

which implies G = H = 0. 

The proof of (j60p is an easy modification of the above argument based on induction, 
duality and density arguments as done in [29]. The tangential regularity can be obtained in 
the same way. The only difference for the normal regularity is that the degeneracy of the 
linear operators L^^C^ in ()58p will change according to the number of d-^ for G and thus the 
appropriate function space is X'^'^ as in the same spirit of Section 15.31 We omit the details. 

5.5 Formal derivation of the equation for G and H = Curl^G 

We first derive the equation for F = w''^ 83(1x1^^°' 837]) = wd'^rj + (1 + 0)831x1837]. From (j35p . 
we obtain 



9?^' + ^ (w'+'' {w8l[A'^J-^/"] + il + a)83w83[A'^J-'/'^]} ,^ ,) ,k 

- (1 + a)w83W,k83[A'^J-^h +wl°'^ + (1 + 0)83x0!^'^ = 

First we look at (**) 

u;9f [4*^J-^/"] + (1 + a)83w83[A^J-^^''] 
= -{J-^/^AlAl + i J-V°4'=^^}{u;a|r?^, +(1 + a)83W 83V' ,s } 

- w83{J-'/''A';At + i J-V-^A?} . ^3^'^,^ 

= -{j-^/"A';At + i J-^/°4'=^,^,}F^, 

+ { J-i/°A^A| + 1 J-i/°Af A?}{u;„ 8lrj' + (1 + a)83W,s ^sr/l 

- w83{J-'/''A',At + 1 J-V-Af A?} . 83r]\s 

Now we rewrite (•) as follows: 
1 



{*) 



1 



W ■ W 



i+"{u;5|[Af J-1/"] + (1 + a)a3u;93[Af J"i/"]} , 



w^ 



w 



l+a 



{^ai[Af J-1/"] + (1 + a)83w83[A^,J-^l''\}) ,k 



+ w,k {w8l[A'^J-'/"] + (1 + a)83W ^si^f J-^/"]} 

+ ^ («^'+"{ J-'/"vl,^4' + i J-^/°Af A?}{^„ afr?^ + (1 + a)93«^,s 537/'-}) ,fc 



+ u;,fc {u;ai[-4f J-^/"] + (1 + a)83w83[AU-'h} ^..^ 



-(*) 
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We now show that the undesirable terms wd^rj and d^rj cancel out in (*) + (**) and thus the 
dominant terms are of wd'^'r],^ and dsri,a-- Indeed, (*) + (**) can be written as follows: 

(*) + (**) 

= W-- (u^+"{J-'/"A^AI + ^J~'/"AfAf.}{w,s dill' + (1 + a)dsw,s dsv'}) ,. 

+ wJ-'/"{[A^,A^ + hAfA'^,]w,, div' - [A^AI + ^AfA'^jw,, dlrj',^ -[A';A^ + \A\AI\ 

■ w,k 9|77^. +[A^,AI + i^^?] [(2 + a)d3W,s dlr,' + (1 + a)dlw,s ^37?^} 
+ (1 + a)«;,3 J-^/'^ilA^AI + ^AUr^.. dlv"- - [A,«A| + ^AtA',]^,^ d^r,\s 

- [A^,A1 + ^A'^r]^,fc a37?^, +(1 + a)[AlAl + ^AlAl]d^w,s ^sr?"} 

+ w{d^[J~^/''AlA1 + ^J-^'^AlA'^,]w,„ dh'' - ds[J~'/''A';At + ^J~'/''A^A',]w,^ d^r^\s 

- d^\J-^AlA1 + 1 J-Uf A-]u;,fc d^rj\„ +(1 + a)d^[J'^ AIAI + ]^J-"^AlAl]d^w,s d^rj'} 

Therefore, the equation for F reads as follows. 

d^tF' - w;-"(u;i+"{J-i/"A^4^ + ^J~^I^A\AI}F\s ),k 

- w-'^iw^+^d-iiJ-^l'^AlAl + i J-V^yi^^A^} . a3r?^, ),, 

+ w-^ (w^+^{J-"^A-,At + iJ-^/^AIADiw.s dh' + (1 + a)dsw,, d^v'}) ,« 

+ wJ-'/^iA^.A'l + iAfA-,]w,^ dir,' - [A^A^ + iAfA',]w,^ dirj' ,s -[A^AI + i^lf A^] 

• w,k dill"-,, +[AlAt + ^AUI] [(2 + a)d:,w,s dh' + (1 + a)dlw„ 83^!} 
+ (1 + a)w,3 J-V"{[^3^- + ^4^::]^/;,, afr?'- - [A'^Af + i^r^?]w^,« 93r?^. (64) 

- K^Af + ^Af^^]u;,fe d3v',a +(1 + a)[A34^ + ^AfA',]d3W,s d^v'} 

+ u;{a3[j~i/"A34- + i j-i/^Af^^ju-,,, afr/'' - a3[j-^/"A;:A| + ^j-^/'^AfAf.] 

■ w,^ dsv'^s -d3[J-^A>;A1 + ^J-7^A^A'^]w,k dsv'',^ 
+ (1 + a)d3[J-iAfAI + lJ-^A^Af.]dsw,s d^vl 

- (1 + a)u' d3W,k dslA'lJ'^^'^] + w /"'^ + (1 + a)d3W 1°'^ = 

The first line is the main part, the second line has full derivative but with the desirable 
weight w'^, the rest of lines are either of lower order with respect to ^3 with appropriate 
weights. We denote the second line through the last line by R^ . Next by using the fact 
A'l = S^Al = A^rf ,s Al, we write the equation (fT2]) as follows: 

Now combining (|64p and (j65p with the equation for [dl + (?|]r/ in ([37]), the equation for 
G = [—dl — i9| — i(;^"93it;^+"53 + A]r/ reads as follows: 

d'^G' - u'~"(u;^+"{ J-^/"^^4" + ^ J-l/"A,^^^,}G^, ),fe 

= -w-^{t,^+'^dp[J~^'^A^,At + i J-i/"^fc^^] . g^^'^^^ )^^ +j2,o,o (66) 

+ ff - A(2 + ^)u;-"('u;i+" J-^/"^f ),fc = 7^^ 
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We note that TZ^ consists of essentially lower-order terms and in particular, 
|l7^*||xa,Jv-2 < oo if £^ iVi^^tV) < oo and ||-C'u'||;5ca,iv < oo . 
The equation for H = Curl^G can be derived in the same way from the curl equation (I42p . 
dtH = [-dj - wdl - (2 + a)d:iw ds + X][Cm\uo]'j 

- ^tA'^,^l + wdl + (2 + o)d^wd^ + A]7?^, -dtAi\dl + wdl + (2 + a)d^wdj, + a]V„ 

p=i 

- [f^Al ■ dtdl^r^^,s -wdilAi . dtdl-^7f,s ) 

+ (2 + a)dMdzA'j ■ dtv\s -dsAi ■ dtv^ ,s ) + A / {dtA'j ■ Otv^s -dtA^ ■ dtrf ,s )dT 

Jo 
2 ,.t (67) 

- E / (^*^^j' ■ ^tdl'^v^s +wdtdlA] ■ dtdl-^7]\, 

- dtd^^Ai ■ dtdl-^rf^s -wdtcf^Al ■ dtdl-^Tj^, )dT 

1 ft 

- V(2 + a)d^w / [dtdlA] ■ dtdl-''^\s -dtdlAl ■ 9*5^^^,, )dr 
p=o -^0 

+ d^wdtiA'^j ■ d^v^a -Al ■ d3rf,a ) - dt{w,^ A] ■ dlri' - w,^ A^ ■ dlrf) 

- (2 + a)dt{d3W,s A'j ■ dsv'' - d3W,s A| • O^Tf) = [S]'^ 

We note that 

||5||2c«,iv-2 < OO if £■ {r],dtri) < oo and ||Dw||j!(^a,iv < oo . 

6 General smooth initial domain Q 

Here we would like to discuss the changes to be made to the argument to prove theorem 12.21 
for general domains. 

There are different ways of trying to extend the result to the general case. One can use 
K charts to cover the boundary of fi and one chart for the interior. Then, one can use 
change of coordinates to straighten out the boundary for each chart. One can then prove a 
priori estimates. However, we think that one of the disadvantage of this method is that the 
proof of existence of approximate solutions is technical since one has to solve K + 1 problems 
simultaneously at each step. Here, we will present a more geometric way motivated by Shatah 
and Zeng |58]. 

Recall the notations introduced in Section [2^ The main difference between the tangential 
and normal derivatives is the fact that |i9fl^w| < Cw and that ^ < d(^w < C. 

The proof of the a priori estimates, namely Proposition 14.11 is identical. One has just to 
replace d^ by 5^ and di , 82 by dp for (3 e T. 

When solving the iteration scheme, we replace the definition of G by 

G = [Y^ {dp)*dp - w-'^d^w^+'^d^ + X\u = Qu (68) 

/3er 
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Jl^\Fi , ..li/^\P^ , ^Pl 



we recall that if dp = a]^{x)di + 03 (x)92 + ^3(^^)53 then (9/3)*(.) = 9i(a]^(x).) + 52(02 (^)-) + 



93(a3(x).). Actually, a better choice of the operator Q is to take it to be selfadjoint in 



L (w°'dx) and hence take 

G=[Y^ w-'^idpfw^'dp - u;-°(a^)*u;^+°a^ + \]u = Qu. (69) 

/3er 

As in section ISTTl we solve the following approximate system for G^+i and H^j^i: 

w^'d^G.+i + jzt G.+i + z:^ G.+i + /:^ i^.+i = w'^n, 

where TZ^, and S^ are given by by formulae that are similar to those in section 15.11 and consist 
of lower order terms. The rest of the proof is identical. We just mention that Lemma 15.21 
was proved above in the flat case. Extending it to the case of the operator Q defined in (|68p 
or ()69p can be easily done following the arguments in [5li6|. 

7 Discussion 

The study of vacuum states in gas and fluid dynamics can be tracked back at least to a con- 
ference, Problems of Cosmical Aerodynamics, held in Paris, 1949 where one session chaired by 
J. von Neumann was devoted to the existence and uniqueness or multiplicity of solutions of 
the aerodynamical equations [60]. See also the review article by Serre [57J. The last question 
raised by von Neumann concerns the validity of the Euler system in the presence of a vacuum: 

There is a further difficulty in the expansion case considered by Burgers. It was accepted 
that the front advances into a vacuum. It is evident that you cannot get the normal condi- 
tions of kinetic theory here either, because the density of the gas goes to zero at the front, 
which means that the mean free path of the molecules will go to infinity. This means that if 
we are in the expanding gas and approach the (theoretical) front, we will necessarily come to 
regions where the mean free path is larger than the distance from the front. In such regions 
one cannot use the hydrodynamical equations. But, as in the case of the shock wave, where 
ordinary conditions are reached at a distance of a few mean free paths from the shock itself, 
so in the case of expansion into a vacuum, at a short distance from the theoretical front, 
one comes into regions where the mean free path is considerably smaller than the distance 
from the front, and where again the classical hydrodynamical equations can be applied. If this 
is applied to expanding interstellar clouds, I think that in order that the classical theory be 
true down to 1/1000 of the density of the clouds, it is necessary that the distance from the 
theoretical front should be of the order of a percent of a parsec. 

This statement led to interesting discussions, which are still enlightening at present, 
among von Neumann and several participants including Heisenberg. An important com- 
ment made by Heisenberg was the role of boundary conditions related to the boundary layer 
theory. And von Neumann answered: 

The boundary layer theory for a fluid of low viscosity certainly furnishes a monumental 
warning. The naive and yet prima facie seemingly reasonable procedure would be to apply 
the ordinary equations of the ideal fluid and then to expect that viscosity will somehow take 
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care of itself in a narrow region along the wall. We have learned that this procedure may lead 
to great errors; a complete theory of the boundary layer may give you completely different 
conditions also for the flow in the bulk of the field. It is possible that the same discipline will 
be necessary for the boundary with a vacuum. 

It is amazing to learn that the difficulty and importance of the problem with vacuum was 
already discussed more than 60 years ago. Although the mathematical theory of vacuum and 
related subjects in the domain of kinetic theory and boundary layer theory is far from being 
complete even now, it is very exciting to make a first step towards a better understanding of 
the problem. 

Our methodology developed in this paper turns out to be robust and it would shed some 
light on other vacuum free boundary problems in more general framework such as problems as 
discussed in Section [L2l We believe that our energy method can be adapted to the relativistic 
Euler equations in a physical vacuum [32]. And in the upcoming paper [31], with this new 
approach, we will give different proofs of the well-posedness of compressible liquid in vacuum 
studied in [38] and the well-posedness of smoother vacuum states of compressible Euler flow 
considered in 1451. 
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